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ABSTRACT 

> ■ 

We present an analysis of weak lensing observations for RXJ 1347-1 145 over a 
. 43' x 43' field taken in B and R filters on the Blanco 4m telescope at CTIO. 

RXJ1347-1145 is a massive cluster at redshift z = 0.45. Using a population of 
galaxies with 20 < R < 26, we detect a weak lensing signal at the p < 0.0005 
level, finding best-fit parameters of o~ v = MOOI^q km s _1 for a singular isothermal 
^p" 1 . sphere model and r 2 oo = 3. 51^2 Mpc with c = 15lfjj for a NFW model in 

an Vt m = 0.3, = 0.7 cosmology. In addition, a mass to light ratio M/L R 
= 90±20 Mq/Ljio was determined. These values are consistent with the previous 
weak lensing study of RXJ1347-1145 by Fischer & Tyson (1997), giving strong 



evidence that systemic bias was not introduced by the relatively small field of 
view in that study. Our best-fit parameter values are also consistent with recent 
X-ray studies by Allen et al. (2002) and Ettori et al. (2001) but are not consistent 
with recent optical velocity dispersion measurements by Cohen & Kneib (2002). 

Subject headings: gravitational lensing — galaxies: clusters: individual (RXJ1347- 
1145) 
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1. Introduction 

In this paper, we present recent weak lensing observations and analysis of RXJ1347- 
1145, the most luminous galaxy cluster in the ROSAT All Sky Survey (Schindler et al. 1995), 
which lies at a redshift of 0.45. RXJ1347-1145 has been the subject of numerous observations, 
including a previous weak lensing study by Fischer & Tyson (1997), and several X-ray and 
optical studies. The primary purpose of the wide-field weak lensing observation presented 
in this paper is to test and confirm the earlier gravitational lensing observations of Fischer 
& Tyson (1997). 

Early gravitational lensing studies of RXJ1347-1145 determined a different total mass 
than early X-ray studies. Such observed disparities for a large and luminous cluster like 
RXJ1347-1145 raised the possibility of potential systemic problems in methodology. As 
galaxy cluster mass distributions are a useful probe for dark matter models, the lensing, 
X-ray, and optical communities independently sought better observations of this cluster. 

Fischer & Tyson (1997) studied RXJ1347-1145 through weak lensing analysis of data 
taken on the Blanco 4m telescope at CTIO using the prime focus CCD camera in 1995. 
Using an approximately 14' x 14' field of view, they found that the mass interior to a 1 
Mpc radius was 1.7 ± 0.4 x 10 15 M Q . Assuming an isotropic velocity distribution, this mass 
corresponded to a velocity dispersion of 1500 ± 160 km s -1 . 

The results of Fischer & Tyson (1997) were not consistent with the X-ray results of 
Schindler et al. (1997), who concluded that the temperature of the cluster was 9.3lJ;o keV. 
As their study presented no evidence that there was a temperature drop from the center of 
the cluster, Schindler et al. (1997) assumed a constant temperature to find that the total 
mass within 1 Mpc was 5.8 x 10 14 M Q . 

There were several, independent concerns regarding these weak lensing and X-ray studies 
that led to further observations. The primary concern for the X-ray studies was the constant 
temperature assumption in the presence of an apparent large cooling flow. For the weak 
lensing studies, there was a possibility that the relatively small field of view would introduce 
a systematic bias, as the outer weak lensing measurements directly influence the predicted 
interior mass. 

Since 1997, several X-ray studies of RXJ1347-1145 have been carried out, most recently 
Chandra (Allen et al. 2002) and BeppoSAX (Ettori et al. 2001) observations. The Bep- 
poSAX observations of Ettori et al. (2001) determined that the 3a lower bound of the X-ray 
temperature of RXJ1347-1145 was 12.1 keV, ruling out the 9.3TJl;o keV temperature measure- 
ment of Schindler et al. (1997) at the 99.99 % confidence level. The Chandra observations of 
Allen et al. (2002) determined the temperature of the cluster to be 12.2 ± 0.6 keV; in rough 
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agreement with the predictions of Ettori et al. (2001). 

This temperature measurement of Allen et al. (2002) is consistent with the analysis 
of ASCA data by Allen & Fabian (1998). The temperature is also confirmed by two mea- 
surements of gas temperature using the SZ effect - those of Pointecouteau et al. (2001) 
and Komatsu et al. (2001). The mass modelling consistent with Allen et al. (2002) is also 
consistent with the weak lensing observations of Fischer & Tyson (1997). 

The recent X-ray measurements of Allen et al. (2002) are not, however, in agreement 
with the recent Keck spectroscopic survey conducted by Cohen & Kneib (2002). Surveying 
47 spectroscopically confirmed members of RXJ1347-1145, Cohen & Kneib (2002) determine 
a central mass that is significantly lower than the X-ray or previous lensing studies. Cohen 
& Kneib (2002) summarize, in some detail, the current X-ray and strong lensing predictions 
and the past weak lensing predictions in Table 4 of their paper. 

The weak lensing observations presented in this paper are based on images of an ap- 
proximately 43' x 43' usable field of view taken in both the R and Bj bands on the Blanco 4m 
telescope at CTIO in 2000 with with light-to-moderate cirrus and 0.8 — 1" seeing. Although 
the imaging depth of this study is less than that of Fischer & Tyson (1997), the observations 
presented here include an approximately nine times larger field of view with better seeing. 

The primary goal of this wide-field gravitational lensing study is to rule out any systemic 
bias introduced by the smaller fields of view in Fischer & Tyson (1997). A secondary goal of 
this paper is to investigate the discrepancy between the current X-ray and optical velocity 
dispersion measurements of RXJ1347-1145. 

Because the observations presented here were designed for weak lensing analysis, they 
do have sufficient angular resolution to definitively judge either the X-ray or optical studies, 
which focus in detail on the central region of the cluster. However, mass models from this 
wide-field weak lensing observation can be compared with other recent studies and will add 
to the overall understanding of RXJ1347-1145. 

2. Observations 

The observations were carried out on March 10, 2000 using the Blanco 4-meter telescope 
at the Cerro Tololo Inter-American Observatory with the MOSAIC II 8k x 8k camera. 
A total of ten 720-second R and eight 420-second Bj exposures were taken with light-to- 
moderate cirrus and 0.8 — 1" seeing. The usable field of view was approximately 43' by 43', 
given the pixel size of 0.258" and a large dithering pattern needed for good night sky flats. 
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The images were processed using the IRAF package mscred up to the flatfielding; the 
image registration and stacking was made with a preliminary version of the Deep Lens 
Survey (DLS) pipeline software (Wittman et al. 2002), which convolves the images with a 
circularizing kernel (Fischer & Tyson 1997) during the combining step, resulting in a circular 
point spread function (PSF). Object catalogs were generated with SExtractor (Bertin & 
Arnouts 1996), and the shape parameters were re-measured using ellipto (Bernstein & 
Jarvis 2002) for all the objects. 

Because the night was not photometric, we were unable to obtain a precise zeropoint. 
Furthermore, we had too few unsaturated USNO stars to permit zeropointing by matching to 
the USNO-A catalog. We were able to achieve a rough calibration by matching photometry 
with that for objects from the previous weak-lensing observations of Fischer & Tyson (1997), 
which used the same telescope (but a different camera). Because in this paper we are 
interested in the shape information primarily, and the photometry and colors are used only to 
construct relative color-magnitude diagrams for object selection, our results are not affected 
by this zeropoint ambiguity. However, we caution against using the data presented here for 
absolute photometry. 

3. PSF Anisotropy 

In general, a spatially variable intrinsic point spread function anisotropy (PSF) will be 
present in our image. This anisotropy must be removed before determining the weak-lensing 
induced shear. We corrected for PSF anisotropy by finding fourth-order polynomials of the 
spatial position that circularize the PSF under convolution with the original image, as in 
Fischer & Tyson (1997). 

To determine the convolution kernel, we measured the initial moments, {I xx , I yy , I xy }, 
for the stars in our image and solved for a kernel that, when convolved with the original 
image, yielded star moments l£ x = I* , and J* = 0. 

Potential stars were selected from objects in the box of Fig. 1. Because the PSF size 
varied vertically across the image, any single R magnitude versus size criteria either left 
regions of the image underpopulated or selected too many galaxies. Therefore, potential 
stars were selected using an R magnitude versus size criterion that varied across the image 
to ensure that the star population contained members from all regions. Note that this 
necessarily implies that many of the objects in the box of Fig. 1 were not included in the 
initial star selection. 

As shown in Fig. 2, this potential star sample was contaminated by a large number of 
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unresolved, bluish galaxies. Therefore, a color cut was made to exclude objects with observed 
color Bj — R > 2.0. After this color selection, we obtained a sample of 1063 non-saturated 
stars, or approximately 0.57 stars per square arc minute. 

The average (uncorrected) ellipticity, 

/- - \ _ / [xx - lyy 21 xy \ /i\ 

(ei, e 2 ) - ( ), (1) 

\ J xx i J yy *xx ~ 1 yy / 

for these stars was 



e 1 = -0.0030 ±0.0008, 

e 2 = 0.0010 ±0.0005. (2) 

This sample had an average PSF size of (1/2) (7 XX ± I yy ) = 0.631 ± 0.003 square arc seconds. 

Applying our spatially dependent, fourth-order convolution kernel to our original R 
image was successful in circularizing the stars in our sample. After the convolution, our 
1063 star sample is characterized by ellipticities 



e 1 = 0.0007 ±0.0008, (3) 

e 2 = 0.0000 ±0.0005. (4) 

The star sample's PSF size increased slightly to 0.639 ± 0.003 square arc seconds, as is 
expected under a circularizing convolution. 



4. Weak lensing analysis 

In this section, we discuss issues related to object selection, systematic errors, and weak 
lensing detection. We consider two object selection criterion, and show that by maximizing 
the number of objects available for analysis, a highly significant weak lensing detection is 
made that is not influenced by systematic errors. 



- 6- 



4.1. Object selection 

Selection of potential galaxies in the background of RXJ1347-1145 was achieved using a 
variety of cuts used to eliminate stars and saturated, poorly resolved, or foreground objects. 
Such objects would reduce the strength of the weak lensing detection because their shape is 
not correlated with the gravitational lensing of the cluster. 

First, the shapes of all objects detected by SExtractor were determined by the program 
ellipto, which is discussed in Bernstein & Jarvis (2002). Ellipto carefully separates objects 
and measures moments based on user defined criteria beginning with the shape and centroid 
measurements from SExtractor. Because the seeing in the Bj image was poorer than the 
R image, the shapes measured by ellipto were based on the R image only. 

Using a set of difficult to deblend objects from our image, we determined that to achieve 
optimal shape measurements, the centroid of a given image could not be allowed to wander 
more than 30% of the size of the object as ellipto determined the object's shape. Objects 
whose shape ellipto could not determine were excluded from further analysis. 

Figure 3 shows a plot of size versus magnitude of the objects whose shape ellipto 
could determine. Saturated objects appear in the arc along the left hand side, while most 
stars appear in the leftward arm. We see in Fig. 4 that almost all of the objects that can be 
distinguished from the stars by a size cut 0.5 x (I xx + I yy ) > 7 have R magnitudes between 
18.0 and 26.0. 

In principle, the color, c = Bj — R, of the cluster can be used to eliminate foreground 
galaxies by selecting only those objects whose color is more blue than the cluster. However, 
preliminary analysis showed that color cuts generally reduced the strength of the weak lensing 
signal in this study. 

The reason that color cuts reduced the measured shear in the current study was two- 
fold. First, because the seeing was relatively poor in the Bj images, approximately 10% of 
the objects in our sample did not have well determined Bj magnitudes, even though they had 
well measured shapes in the R image. Hence, any color cut excluded a substantial number 
of objects in the background of the cluster whose shapes were well measured. Second, while 
color cuts do eliminate most foreground objects, they also eliminate some background objects 
that contribute to the shear. 

Since we can not rely on color selection, the primary object selection criteria of im- 
portance will be the R magnitude cut. We present in this section and in our analysis two 
different R selections. 

It would appear that the conservative way to proceed would be to make an R magnitude 
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cut above the "expected" R magnitude of cluster members at redshift 0.45. This would seem 
to ensure that potential cluster alignment does not contaminate the weak lensing signal. 
Based on data from Fukugita et al. (1995), cluster members at z — 0.45 have expected 
fT-corrected Bj magnitudes of 23, and the expected Bj — R color of such members is 
approximately 2.7. (This agrees well with a sample of potential cluster members hand 
selected from our image.) Therefore, an R magnitude cut of R > 22 conservatively selects 
objects likely to fall behind the cluster. 

The more aggressive way to proceed is to choose a lower R magnitude cut that maximizes 
the total detected gravitational shear. This is reasonable because, while the R > 22 cut does 
exclude most objects in the foreground of the cluster, it also excludes objects behind the 
cluster that do contribute to the gravitational shear. 

In Section 5, we estimate the redshift of each potential background galaxy using statistics 
from the California Institute of Technology Faint Galaxy Redshift Survey (Cohen et al. 2000). 
Applying statistics drawn from that survey to our sample, we expect that 78.7% of objects 
whose R magnitudes lie between 20 < R < 22 fall behind the cluster, and overall, 90.1 % of 
all objects with 20 < R < 26 are expected to fall behind the cluster. For the more restrictive 
R magnitude cut at R > 22, the overall number of objects expected to fall behind the cluster 
is 93.6 %, a slight increase from 90.1 % for the 20 < R < 26 cut. However, the 22 < R < 26 
cut excludes 22.2 % of all objects with 20 < R < 26 expected to fall behind the cluster, 
significantly increasing the noise. 

Our preliminary analysis determined that a simple R magnitude cut at R > 20 was more 
effective at maintaining a high number of background objects than color cuts or cuts based on 
the expected cluster luminosity function. These findings are consistent with the preliminary 
analysis of Khiabanian (2003), who has analyzed DLS deep field data for combinations of 
parameters that lead to the most robust weak-lensing detections. 

In the analysis sections of this paper, we present best fit model parameters for both the 
20 < R < 26 and 22 < R < 26 cuts with a size restriction 0.5 x (I xx + I yy ) greater than 7. 
However, we feel more confident in the 20 < R < 26 cut as shown in the U-shaped region 
of Fig. 3. Using this constraint, we obtained a well- measured background galaxy sample of 
19796 galaxies, or approximately 11 galaxies per square arc minute. 

4.2. Systematic errors and weak lensing detection 

Residual PSF anisotropy or inherent background galaxy ellipticity correlations (indepen- 
dent of any lensing) would lead to systematic errors in the measurement of the gravitational 
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shear. In this subsection, we present evidence that any remaining PSF anisotropy or inher- 
ent background galaxy ellipticity correlations are minimal in strength compared to the weak 
lensing signal. 

The combined R band image of RXJ1347-1145 used in this study is shown in Fig. 6. 
The cluster is located slightly northeast of center in the frame. Figure 7 shows the two- 
dimensional mass map determined by the ellipticity correlations in the 19796 potential back- 
ground galaxies for 20 < R < 26. This mass map has the same scale as the R band image, 
and was produced by f iatmap, a program that determines a two-dimensional mass profile 
from detected correlations in object ellipticity (Wittman et al. 2002; Fischer & Tyson 1997). 

The color contrast mass map shows a strong (bright) weak lensing signal at the location 
of the cluster. The field appears relatively empty aside from a concentration of mass at the 
location of the cluster. 

Using the center of the cluster determined from the mass map, we determine the average 
gravitational shear in radial annuli centered on the cluster. This task is accomplished by 
determining, for each object in a given radial annulus, the alignment of the object to a circle 
centered on the cluster and then averaging these values. For each annulus, we obtain an 
average tangential shear, j t , (a measure of how aligned objects are to the circle) and an 
estimated error in the tangential shear. 

Figure 8 shows the tangential shear detected for RXJ1347-1145 with la error bars from 
analysis of the 19796 background galaxies with 20 < R < 26. Also shown is the contribution 
of the residual PSF anisotropy to the tangential shear as determined by the stars. We see 
that the residual PSF anisotropy is consistent with zero shear detection and substantially 
weaker than the detection based on the potential background galaxies. 

In addition to measuring average tangential shear, we can measure an average shear 
oriented at 45° from tangent to the circle centered on the cluster. In weak lensing, this 
45° shear should be consistent with zero shear detection at the large radii of this study. In 
dashed lines, Fig. 8 shows the 45° shear detected from the 19796 background galaxies with 
20 < R < 26. That the 45° shear is consistent with zero gives another indication that major 
systematic errors have been controlled. 

As a third test of systematic error, we show the results of an x — y scramble test in 
Fig. 9. In this test, we consider the quintuplet (x\ y\ P xx , P yy) P ) for the ith object used in 
the shear detection of Fig. 8 and decouple the position (x l ,y l ) from the shape information 
{P xx , Iyy, P xy ) by randomly assigning shape information to positions, forming the quintuplets 
(x\ y\ I xx , P , P x ), where r represents the same random integer for each moment. 
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Since weak lensing produces a position dependent shearing, the quintuplet (x\ y\ P xx , I yy , I xy ) 
should result in tangential shearing, but the scrambled quintuplet (x l , y\ I xx , I yy , I xy ) should 
not. In Fig. 9, we see that the x — y scrambled tangential shear is consistent with zero signal. 

We statistically test the null hypothesis (that the true value of 7 t is zero) by computing 
a x 2 statistic assuming that the error in the true value of 7 t is equal to the error in the 
measured This results in a x 2 °f 42.0 for 15 degrees of freedom for the sample with 
20 < R < 26. The probability of obtaining a x 2 this large in the absence of a weak lensing 
effect is less than p = 0.0002. The more "conservative" R magnitude cut at 22 < R < 26 
yields x 2 — 28.5 for 15 d.o.f., which corresponds to a probability of null detection less than 
p = 0.02. 

For comparison, a null hypothesis test for the x 2 statistic for the 45° shear shown in 
Fig. 8 results in a x 2 of 16.6 for 15 degrees of freedom, while the tangential shear measurement 
of the stars had a x 2 °f 26.6. The x — y scramble test yields a x 2 of 11.4 for 15 degrees of 
freedom, and is consistent with a null detection at p — 0.724. 



5. Mass reconstruction 

In gravitational lensing, the tangential shear, 7 t , at a given projected radius r in the 
lens plane is related to the surface mass density, E, at that radius. Defining k = E/E cr j t , the 
ratio of the surface mass density to the critical surface density for multiple lensing, we have 



lt{r) = «(< r) - «(r), (5) 

where k(< r) is the normalized, average surface mass density interior to r and R{r) is the 
mean density at r (Miralda-Escude 1991). For a given mass model that determines E, Eq. 5 
allows us to relate the model parameters to the observed shear, j t . 

The critical surface density is given by (Schneider et al. 1992) 

= (6) 
cnt AnG DiDi s 1 ; 

where D s , D\ and D\ s are the angular-diameter distances between the observer and source, 
observer and lens, and lens and source, respectively. Hence, E crit is both redshift and cos- 
mology dependent. 



We will determine an average E cr j t for our sample by using the known redshift of the 
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cluster, z = 0.4509 ±0.003 (Cohen & Kneib 2002), and an average value of T, crit for each the 
many sources. 

To determine the redshifts of the sources, we develop a redshift to R magnitude relation 
using the California Institute of Technology Faint Galaxy Redshift Survey (Cohen et al. 
2000). Dividing the well measured objects in this survey into 0.50- wide R magnitude bins 
from 19.5 < R < 24.5, we determined an initial average value of z for a given R magnitude. 
To remove outliers from each bin, we excluded objects whose measured redshift was more 
than one standard deviation from the bin average, and then recomputed a final average 
redshift for the given R magnitude. 

Plotting this average redshift of each bin against R magnitude, and using a weighted 
least squares method, we determine that R is related to z in a linear fashion 

z = m x R + b, (7) 
where the parameters m and b are given by 

m = 0.094 ±0.055, 
b = -1.5 ±1.2. (8) 

The slope- intercept parameter space is degenerate, so the error bars on m and b from Eq. 8 are 
linked. Figure 10 shows the binned data from the CalTech Faint Galaxy Redshift Survey, 
the best fit line from Eq. 7 and two bounding curves corresponding to a 68% confidence 
interval. 

The Faint Galaxy Redshift Survey does not contain a sufficient number of galaxies 
fainter than R > 24.5 to rule out a possible selection bias. Therefore, we posit that the 
linear fit found from the 19.5 < R < 24.5 region of the survey will apply to our entire galaxy 
sample in Fig. 3, even though our sample contains objects fainter than R > 24.5. 

We note that any error introduced in the redshift by our linear fit for objects with R 
magnitude greater than 24.5 will have a minimal impact on the eventual determination of 
S crit for two reasons. First, from Fig. 4, we see that relatively few objects will be used at 
R > 24.5. Second, almost all objects this faint will lie in the vicinity of z ~ 1 where the 
angular diameter distance flattens as a function of z. Hence, even relatively large errors in 
z around z ~ 1 do not propagate strongly into errors in D s , Di s or S crit . 

Applying Eq. 7 to our sample with 20 < R < 26, we determine that our sample has 
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an average redshift, z = 0.74 ± 0.13, and that the average critical surface density is £ cr j t = 
82OOl^/iM pc~ 2 for an Vt m = 1, fi A = cosmology with # = 50 km s' 1 Mpc" 1 . (Here, 
the upper and lower error limits are estimated using the bounding curves from Fig. 10.) For 
an Vt m = 0.3, Vt^ = 0.7 cosmology, the average S crii is given by 6000+230° hM & pc" 2 . For the 
sample with 22 < R < 26, the values for S crit are 8OOOl^°/iM pc" 2 and 65OOl 2 £° ) ° ) /iM 
pc~ 2 for Q m = 1, Q\ = and f2 m = 0.3, = 0.7 cosmologies, respectively. 

In principle, the critical density varies as a function of angular radius because the 
magnification by the lens of background sources makes sources closer to the center of the 
lens appear relatively brighter than sources at a greater angular radius. This increase in 
brightness causes us to underestimate the redshift of the source and overestimate the critical 
density. 

Fischer & Tyson (1997) used Monte Carlo simulations of background galaxies magnified 
by an isothermal lens with a velocity dispersion of 1460 km s" 1 to determine a plot of 
£ crit versus radius. Figure 7 of their paper shows that the observed shear of their study's 
innermost data points were affected by the variation in S crit . However, accounting for this 
variation did not lead to a significant change in model fit, primarily because the innermost 
angular bins have the fewest data points and hence the highest uncertainty. 

In the current wide-field study, the innermost radial bin falls over the range in which 
Fig. 7 of Fischer & Tyson (1997) shows a significant variation in S crii . However, the remaining 
radial bins fall into a region where the variation of £ crit is within the error introduced by 
using the fit of Eq. 7. 

Therefore, in this wide-field study, we will use a constant value of S crit as a function 
of angular radius for a given cosmological model. We feel justified in doing so because the 
current wide-field study covers almost twice the radial size of Fischer & Tyson (1997) and 
hence is less susceptible to any biasing introduced by errors from the innermost data point. 
These errors are small because the weighted x 2 algorithms used to fit mass models discount 
data points with large errors. 

6. Mass models 

In this section, we apply two mass models to our observed shear of RXJ1347-1145: a 
singular isothermal sphere model (SIS) and the Navarro, Frenk, and White (NFW) model 
(Navarro et al. 1997). Throughout this section, we assume that the Hubble constant has 
a value of H = 50 km s~ x Mpc^ 1 . For each model, we consider an Q m = 0.3, Q\ = 0.7 
cosmology, and briefly cite the best fit parameters for the Q m = 1, = cosmology in 



order to compare with previous studies of this cluster. In 6.3, we summarize the cluster 
masses inferred from the two models. 



6.1. 



SIS model 



In the SIS model, the surface density is parameterized by the velocity dispersion, a v : 



Here, r is the projected radius in meters in the lens plane from the center of the cluster. 
In the SIS model, S(r) = 2£(r), so by Eq. 5, the tangential shear is related to the velocity 
dispersion by 



Using a standard linear regression fitting algorithm that minimizes \ 2 ■, we determine that 
the velocity dispersion predicted by our observed shear in an Q m = 0.3, Q\ = 0.7 cosmology 
is a v = 1400 km s _1 for the 20 < R < 26 magnitude cut (minimum \ 2 — 14.8 for 14 d.o.f.), 
and 1550 km s _1 for the 22 < R < 26 magnitude cut (minimum \ 2 — 12.3 for 14 d.o.f.) 

Two independent factors contribute to error in the inferred velocity dispersion. First, 
the process of minimizing x 2 carries with it a random statistical error that can be estimated 
using the A% 2 argument of Press et al. (1995). We estimate that these statistical errors 
contribute to uncertainty in the velocity dispersion of a v = lAOOt\lo km s _1 for a 68% 
confidence interval with the 20 < R < 26 cut and a v = 1550^ km s_1 for the 22 < R < 26 
cut. Independently, the errors bars on the lensing critical density, contribute to uncertainty 
in a v as a v = 14001^ km s_1 for 20 < R < 26 and 1550^ km s" 1 for 22 < R < 26. 

Assuming the velocity dispersion of a v = 1400^4q km s -1 , the SIS model corresponds 
to a total integrated mass of 1.4 ± 0.3 x 10 15 M within a radius of 1.0 Mpc. 

For an Vt m = 1.0, Q\ — cosmology, the best fit velocity dispersions are 1500^^q km 
s- 1 (statistical) and 1500^0 km s" 1 (error in critical density) for 20 < R < 26 and 1580±2i° 
km s- 1 (statistical) and 1580^0° km (critical density) for 22 < R < 26. 



lt(r) 



k{< r) — k{t) 




(10) 



6.2. NFW model 



The NFW model depends on a concentration parameter c, and a scale radius r s . The 
scale radius is related to the virial radius r2oo = cr s , which is the radius inside which the 
mass density of the halo is equal to 200p c , where p c = 3 ^J^ is the critical density of the 
universe at redshift z. 

From Eq. 5, we have that 



H NFW {< r) - T.NFwi^r) 
it\r) = = , (11) 

and Wright & Brainerd (2000) derive the relevant expressions for T>nfw{< r ) an d £/v.Fw( r )- 
Because our observational data is specified in terms of annular bins with inner and outer 
radii and r a , we must derive expressions for the average density interior to the annulus 
and the average density in the annulus. The average density in the annulus is given by 



E(r) = E(r ) 



-nn) 



2 ' 



(12) 



while the average density interior to the annulus is given by 



S(r) = S(r D )^^ + nrd^t~ 2 - (13) 
r l Q + rf r A Q + rf 

Hence, our observed tangential shear, 7 t , for a given annulus of inner radius r\ and outer 
radius r Q is given by 




We use Eq. 14 with the expressions for S from Wright & Brainerd (2000) to determine 
the best fit parameters of the NFW model. As the NFW model is highly non-linear, we use 
the Press et al. (1995) implementation of the Levenberg-Marquardt method to fit for c and 
r s , minimizing x 2 - 

For an VL m = 0.3, Vt\ = 0.7 cosmology, we find that c = 15, and r s = 0.23 Mpc, 
or r 20 o = 3.5 Mpc for the 20 < R < 26 cut with a minimum % 2 = 15.3 for thirteen 
d.o.f. For 22 < R < 26, we obtain c = 6.8, and r s = 0.69 Mpc, or r 200 = 4.7 Mpc with 
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X 2 = 10.2. Figure 11 shows the best fit NFW model with the observed gravitational shear 
for 20 < R < 26. 

As for the SIS models, we estimate the statistical errors in the x 2 minimization and 
the error propagated from error in S crit . In the case of the NFW models, the error in £ crit 
introduces a small error in the concentration parameter, c = 15l 2 for 20 < R < 26 and 
c = 6.8^q| for 22 < R < 26 and less than 1% variation in r s for both cuts. 

Figure 12 shows the 68 % confidence region in the c-r 20 o plane for the 20 < R < 26 cut, 
yielding error bars r 200 = 3.5+q; 2 Mpc with c = 15^^ f° r a NFW model in an f2 m = 0.3, 
Ha = 0.7 cosmology. The concentration parameter, c is not well constrained. (The confidence 
region for 22 < R < 26 is similar.) 

We estimate error bars on the total integrated mass using the upper left and lower 
right projections of the confidence interval onto the c-r 2 oo axes. The lack of constraint on 
c yields large error bars on the total integrated mass: 2.7t 2 1 ( i x 1O 15 M within r 2 oo and 
1.41^ x 1O 15 M within 1.0 Mpc in the Q m = 0.3, Q A = 0.7 cosmology for the 20 < R < 26 
cut. 

The best fit NFW parameters for an il m = 1, = cosmology are c = 11.5l£;jj (error 
from Sent), r s = 0.193 Mpc, or r 20 o = 2.22 Mpc with \ 2 — 15.3 for thirteen d.o.f. for 
20 < R < 26 and c = 4.8±g;| (error from S crii ), r s = 0.584 Mpc, or r 200 = 2.78 Mpc with 
X 2 = 1.2 for thirteen d.o.f. for 22 < R < 26. Again, the concentration parameter, c is not 
well constrained by the statistical fit, as is shown in Fig. 12. 

6.3. Observed Mass 

Table 1 summarizes the inferred mass interior to 1 Mpc of the cluster in SIS and NFW 
models for two cosmo logical models for the 20 < R < 26 cut, while Table 2 presents the 
same information for the 22 < R < 26 cut. We see that within a cosmological model, the 
NFW and SIS inferred masses generally agree within la error bars, and that the choice 
of cosmological model makes a small difference in inferred mass. We also note that the 
22 < R < 26 cut generally yields a higher total integrated mass than the 20 < R < 26 cut. 

Figure 13 plots an observed measure of the mass densitometry of RXJ 1347-1 145 com- 
puted by the relation 



5/t(rj) = n{r < Ti) — k{ti <r <r Q ), 



(15) 
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introduced by Fahlman et al. (1994), where r a is the outer radius of the outermost annulus. 
Densitometry measures from SIS and NFW models in an Q m = 0.3, = 0.7 cosmology are 
shown assuming model parameters from best-fits to the gravitational shear data. The error 
bars are statistical error bars computed from the number of galaxies used in the calculation 
only. 

Because the observed densitometry measure uses all of the objects between and r 
to determine 8k at r,^ the values of the observed densitometry measure in Fig. 13 are not 
independent. In fact, small changes in values at the outer annuli give rise to large changes 
in the densitometry value at inner annuli, so that the actual uncertainty in the innermost 
annulus is quite large. Thus, although neither the SIS or NFW model fit the actual value 
of the innermost point, both models are adequate fits to the total profile and yield similar 
total integrated masses. 



We present a rough estimate of the light distribution of the cluster drawn from the R 
band image in Fig. 14. This estimate is limited by the difficulty of accurate background 
subtraction due to the lack of a high quality Bj image that can be used for color selection. 
This luminosity estimate is also limited by the compact nature of RXJ1347 and the relatively 
low angular resolution of the central region in this study. In this section, all measured values 
are the SExtractor values for the combined R band image. 

We begin by selecting all non-saturated objects detected in the R band image by 
SExtractor with a minimal lower limit size cut 0.5 x (I xx + I yy ) > 2 and deleting spu- 
riously detected objects in regions near stars. Because we can not make a color or redshift 
selection, we assumed that all objects are at the cluster redshift, and determined a total flux, 
/, in radial annuli centered on the cluster. From the total flux in each annulus, we determined 
a total absolute magnitude M, applying a K correction of 0.625 from Fukugita et al. (1995). 
A luminosity relative to the solar luminosity in each annulus was then computed from this 
total magnitude. 

To compare this light distribution with the mass distribution determined by the weak 
lensing, we first determined that the luminosity per unit area is well fit by a general SIS 
profile. Since this is the case, we can assume that the luminosity per unit area is linearly 
related to the projected SIS mass density, or that 



7. Luminosity 




(16) 
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where the SIS model for mass density £ is given by the best fit parameters found from the 
weak lensing, or a v = 1400^140 ^ m s_1 f° r the = 0.3, = 0.7 cosmology. We compute 
S(r) in units of solar masses per unit area and find the parameters a and b by minimization 
of 

, / I fl /A-( fl + txS(f)) \ 2 

with 

f n, \ ba v 5a v 

Minimizing x 2 from Eq. 17 has the advantage of finding both the approximate "back- 
ground" luminosity and mass to light ratio simultaneously. Formally, the denominator of 
X 2 should also include a term for uncertainty in L R /A. We neglect this uncertainty because 
we are unable to determine its value without information about the background redshift 
distribution, although it is reasonable that it is approximately constant for all radial bins. 

The best fit parameters, a and b, for seven evenly spaced annuli out to approximately 
2.5 Mpc are a = 1.9 L R& pc~ 2 and b = 0.011 L RQ /M G with x 2 — 10 for five degrees of 
freedom. This implies that M/L R = 90 M Q /L RQ with H = 50 km s _1 Mpc -1 . Figure 14 
plots the luminosity per unit area and best fit scaled mass density (from Eq. 16) along with 
la scaled mass density curves. 

Since we do not have error estimates on the luminosity data, we determine only the 
statistical uncertainty present in the fit using the Ax 2 argument of Press et al. (1995). This 
yields an error ellipse whose projections indicate that the 68% confidence interval for the 
mass to light ratio is M/L R = 91 ± 20 M Q /L RQ . 



8. Discussion 

In this section, we discuss the two R magnitude cuts and compare our results with the 
previous weak lensing results of Fischer & Tyson (1997), the recent X-ray studies of Allen 
et al. (2002), and the recent optical studies of Cohen & Kneib (2002). Our results indicate 
that there were not systematic problems in the previous weak lensing studies. 
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8.1. Object selection criteria 

Because our Bj imaging was poor, color selection criteria dramatically reduced the 
significance of our weak lensing detection by eliminating objects with well measured shapes 
that were undetected in the Bj image. Thus, we were forced to employ only size and R 
magnitude cuts to restrict our image. A size cut (0.5 x {I xx + I yy } > 7) was used to 
eliminate poorly resolved objects and stars, while the R magnitude cut essentially was used 
to control foreground objects (see Fig. 3). 

In principle, it is very important to eliminate foreground and cluster objects from the 
catalog before analyzing gravitationally induced shears for two reasons. First, foreground 
objects weaken the signal because their shapes are not correlated with the weak lensing 
shear. Second, one might conjecture that cluster members align in some correlated way and 
unduly influence the gravitational shear detection if they are included. 

In practice, we find that "restrictive" object selection criteria do not enhance the weak 
lensing detection. In this study, the conservative 22 < R < 26 yielded a substantially lower 
total x 2 relative to the null detection (x 2 = 28.5, 15 d.o.f.) than the 20 < R < 26 cut 
(X 2 = 42.0, 15 d.o.f.). 

The R magnitude cut at R > 22 is reasonable based on the expected cluster R magnitude 
of R 20 at z = 0.45. This cut successfully eliminates most foreground and cluster members 
from further analysis. This would apparently lead to a strong detection of gravitational shear 
by reducing the noise. 

However, a cut at R > 22 also eliminates a significant number of galaxies in the back- 
ground of the cluster that are gravitationally lensed. The effect of removing these lensed 
objects from analysis weakens the signal much more than the signal is enhanced by remov- 
ing the noisy foreground, and so the weak lensing detection is more significant with a less 
restrictive cut at the expected cluster R magnitude. 

We do not believe that inclusion of objects with 20 < R < 22 significantly contaminates 
the sample or biases the weak lensing detection. In Fig. 5, we see that objects in the 
innermost radial bins are slightly brighter than in the outer bins. However, if we consider 
the over-density of objects with 20 < R < 22 in the "cluster" region interior to an angular 
radius of r < 400", we estimate that less than 100 additional objects have been added to the 
catalog of over 1300 total objects with r < 400". 

Further, we observe that the mass inferred from the best fit parameters of SIS and NFW 
models for the sample of objects with 22 < R < 26 is only slightly higher than that from 
objects with 20 < R < 26. In fact, the statistical error bars have significant overlap for the 
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two R magnitude cuts. 

We believe that the results from the 20 < R < 26 magnitude cut represent a better 
measurement of weak gravitational lensing than the 22 < R < 26 magnitude cuts because 
the significance of detection is much higher. However, the best fit model parameters are 
similar, and the basic conclusions of the paper are unchanged for either R magnitude cut. 

8.2. Previous weak lensing observation 

The general results of this paper are in agreement with the results of Fischer & Tyson 
(1997), who report that the total integrated mass within r < 1 Mpc is 1.7±0.4 x 1O 15 M , and 
that the velocity dispersion is a v = 1500 ± 160 km s _1 for SIS model assuming an isotropic 
velocity dispersion and an Q m = 1 cosmology. These values are consistent with our results 
for the same model. 

As is seen in Fig. 9, the wide-field weak lensing observations presented here confirm the 
detection of Fischer & Tyson (1997), which was limited to an approximately 400" radius. 
The current wide-field study continues to detect significant gravitational shear at r > 400"; 
however, the impact this shear has on the model parameters in SIS and NFW models is 
minimal. 

The mass to light ratio of the cluster measured in this study agrees relatively well 
with Fischer & Tyson (1997), who find values of M/L R = 75 ± 20 M Q /L RQ . This is quite 
remarkable considering the difference in methods and difficulty of background subtraction 
in the current study. 

That the current weak lensing study obtains results in agreement with the results of 
Fischer & Tyson (1997) gives a strong indication that the concerns regarding the use of 
smaller images in determining properties from lensing studies were unfounded. The large 
field of view of this study and extremely robust weak lensing detection strongly confirms the 
results of Fischer & Tyson (1997). 

8.3. Recent X-ray observations 

Our results are also in good agreement with the recent Chandra X-ray observations of 
Allen et al. (2002). Using an NFW model, Allen et al. (2002) determine a best-fit scale 
radius of r s = 0.40^^2 Mpc with a concentration parameter of c = 5.87~^\'^. These values 
correspond to a virial radius r 2 oo = 2. 351^33 Mpc and an integrated mass within the virial 
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radius of M 200 = 2.291^ x 10 15 M Q . 

The current weak lensing observations find a similar virial radius, r 2 oo = 2.22 Mpc, and 
a consistent mass M 20 o = 2.3^^ x 1O 15 M for the same cosmological model. Figure 12 
shows the 68 % confidence region for the best fit parameters r s and c for the current weak 
lensing study with 20 < R < 26 and the location of the parameters found by Allen et al. 
(2002). The recent Chandra results are also consistent with a 22 < R < 26 magnitude cut. 

Allen et al. (2002) report that a SIS model does not provide a good fit to the Chandra 
observations. Our observations are more tolerant of a SIS model because we do not have 
lensing observations close to the center of the cluster, a region that is well sampled and 
integral for the X-ray study's conclusions. 

8.4. Recent spectroscopic observations 

This weak lensing study is not consistent with the results of the most recent optical 
velocity dispersion study by Cohen & Kneib (2002). Cohen & Kneib (2002) report a velocity 
dispersion of a v = 910 ± 130 km s _1 , corresponding to an integrated mass of 6.0 x 1O 14 M 
within 1.0 Mpc. This predicted mass is approximately 38% of the mass estimated by the 
current weak lensing observations, and is inconsistent with our results at the 5a level. 

Cohen & Kneib (2002) suggest that their results may be artificially low if RXJ1347-1145 
consists of two clusters caught in the act of merging in a direction perpendicular to the line 
of sight. This flow would have a minimal effect on the wide-field lensing signal analyzed in 
this study because the wide-field imaging does not have sufficient resolution to examine the 
interior region of the cluster. 

9. Conclusions 

In conclusion, the findings of this paper confirm the earlier weak lensing study of Fischer 
& Tyson (1997) and are consistent with the recent X- ray observations of Allen et al. (2002) 
and Ettori et al. (2001). From these findings, we conclude that there were no substantial 
systematic errors introduced in the older lensing study of this system by Fischer & Tyson 
(1997) through the relatively small field of view. Confirming the older weak lensing obser- 
vations of Fischer & Tyson (1997) using a wide- field observation gives substantial evidence 
that weak lensing may not be significantly biased by relatively small fields of view. 

This study's results are in substantial disagreement with the direct optical velocity 
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dispersion measurements of Cohen & Kneib (2002). As the X-ray and lensing results now 
tend to cluster around a velocity dispersion of approximately 1500 km s _1 , it is unlikely 
that the 910 km s _1 measurement of Cohen and Kneib is correct. Cohen & Kneib (2002), 
recognizing this difference suggest that a possible merger within the cluster perpendicular 
to the line of sight could account for the observed difference. 

In the current study, we did not attempt to measure weak lensing signals interior to the 
two arc candidates located approximately 35" from the center of the cluster. Because the 
original goal of this study was to test the small field of view of Fischer & Tyson (1997), the 
wide-field image used in this study does not have enough resolution or depth to conclusively 
measure weak lensing signals in the central region of the RXJ 1347- 1145. 

For this reason, we can not attempt to locate substructure within the cluster to test 
Cohen and Kneib's assertion in this study or compare our results to recent strong lensing 
observations. In addition, we suspect that it is the X-ray observation's greater sensitivity to 
the central region of RXJ1347-1145 that accounts for the inability of Allen et al. (2002) to 
find a good SIS fit. 

Future high resolution lensing studies of the interior region of RXJ1347-1145 may yield 
interesting results. We note that the two-dimensional mass map presented in Fig. 7 does 
suggest some structure, leading to the hope that high resolution weak lensing studies of the 
central region may shed light on possibility of a significant merging event. 

This paper is based on observations made at the Cerro Tololo Intra- American Observa- 
tory which is operated by the Association of Universities for Research in Astronomy, under a 
cooperative agreement with the National Science Foundation as part of the National Optical 
Astronomy Observatories. This work partly supported by NSF grant AST-0 134743. 
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Model Cosmology Parameters Mass (< 1 Mpc) 



SIS 




= 1.0, n A 


= 0.0 


a v 


= 1500T^° km s- 1 


1.6 ± 0.3 


SIS 




= 0.3, Q A 


= 0.7 


a v 


= 14001^ km s- 1 


1.4 ± 0.3 


NFW 




= 1.0, fi A 


= 0.0 


c = 
r s 


= 11.5 

= 1.93 Mpc 


i c+0.6 
i -°-0.5 


NFW 




= 0.3, Q A 


= 0.7 


c = 


= 15.4 

= 2.29 Mpc 


1 4+0-5 



Table 1: SIS and NFW model best fit parameters and total integrated masses within 1.0 Mpc 
(in units of 1O 15 M ) for two different cosmological models with h = 0.5 for gravitational shear 
detected from objects with 20 < R < 26. The VL m = 1 models are consistent with a previous 
weak lensing study and the most recent Chandra X-ray observations. The uncertainties in 
mass assume S cr j t = 8200 h M Q pc~ 2 for the Vt m = 1 cosmology and Tj crit = 6000 hM & for 
the Q m = 0.3, Q\ = 0.7 cosmology. 
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Model Cosmology Parameters Mass (< 1 Mpc) 



SIS 




= 1.0, n A 


= 0.0 


a v 


= 1580^° km s- 1 


1.8 ± 0.4 


SIS 




= 0.3, Q A 


= 0.7 


a v 


= 1550+^° km s- 1 


1.8 ± 0.4 


NFW 




= 1.0, tt A 


= 0.0 


c = 
r s 


= 4.8 

= 0.584 Mpc 


9 9+O.5 
z - z -0.5 


NFW 




= 0.3, Q A 


= 0.7 


c = 


= 6.8 

= 0.691 Mpc 


9 1 +0.2 
0.3 



Table 2: SIS and NFW model best fit parameters and total integrated masses within 1.0 Mpc 
(in units of 1O 15 M ) for two different cosmological models with h = 0.5 for gravitational 
shear detected from objects with 22 < R < 26. The uncertainties in mass assume S crit = 
8000 h M & pc -2 for the Vt m = 1 cosmology and S crit = 6500 h M© for the Vl m = 0.3, fi A = 0.7 
cosmology. 
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18 20 22 

R magnitude 

Fig. 1. — A plot of size versus i?-magnitude for well-defined objects in our sample. To 
determine a convolution kernel that circularized the PSF, potential stars were selected from 
the boxed region using a spatially variable criterion. Note that many objects in the boxed 
region were not included because they do not meet the additional spatial location criteria. 
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Fig. 2. — A histogram indicating the number of objects of a given color in the potential star 
sample in the field of view of RXJ 1347- 1145. The bi- modal distribution indicates that the 
potential star sample might be contaminated by blue galaxies. A color cut of B — R > 2.0 
was applied to limit the sample to likely stars. 
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Fig. 3. — A plot of size versus -R-magnitude for well-defined objects in our sample. Potential 
background galaxies were selected within the "U" shaped region shown for detection of a 
weak lensing signal. 
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Fig. 4. — Two histograms indicating the number of objects of a given R magnitude. The 
upper histogram shows all objects from Fig. 3, while the lower histogram shows only those 
objects whose size, 0.5 x (I xx + I yy ), is greater than 7.0. 
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Fig. 5. — Three histograms indicating the total number of objects (solid), the number of 
objects with 20 < R < 22 (long dashed), and ratio (short dashed) in radial bins. That the 
ratio of the numbers of bright objects and all objects is relatively constant indicates that the 
inclusion of objects with 20 < R < 22 does not significantly contaminate the sample. 
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Fig. 6. — The combined R band image of field around RXJ1347-1145 used in this study. 
The circle, centered on the cluster at an RA of approximately 13:47.30 and Dec of —11:45, 
indicates the region of the image analyzed for weak lensing signals. 
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Fig. 7. — A projected two-dimensional mass map for RXJ1347-1145 over a 43' x 43' field of 
view. The bright, compact, central object is the cluster. The circle indicates the region of 
the image analyzed for weak lensing signals. The relative lack of other substantial regions 
of mass in the middle of the image indicates a significant weak lensing detection about the 
center of RXJ1347-1145. The features detected near the edges of the image are artifacts of 
a reduced amount of data at the boundaries. 
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Fig. 8. — The tangential shear as a function of angular radius due to weak lensing of back- 
ground galaxies near RXJ1347-1145 with la error bars. Also shown is an estimate of the 
contribution to the detected shear due to residual PSF anisotropy as determined by analysis 
of stars (small solid errors bars) and the 45° shear as a function radius (large dashed error 
bars). The Mpc scale is the radial distance in the lens plane for an Vl m = 0.3, il\ = 0.7 
cosmology. 
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Fig. 9. — The tangential shear as a function of angular radius with la error bars is shown 
with the results of an x — y scramble test (dashed error bars). The x — y scramble test is 
a measurement of systemic errors. Also shown are the data from an earlier weak lensing 
detection (triangles). The Mpc scale is the radial distance in the lens plane for an Q m = 0.3, 
f^A = 0.7 cosmology. 
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Fig. 10. — Redshift versus R magnitude data from the CalTech Faint Galaxy Redshift 
Survey. Objects whose redshift fell one standard deviation from an initial bin average have 
been excluded from this plot. The error bars are 1 a error bars for the recalculated bin 
average. The solid line corresponds to a best fit redshift to magnitude relation that is used 
to determine Ti crit . The two dashed lines correspond to 68% confidence limits. 
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Fig. 11. — The tangential shear as a function of angular radius with la error bars is shown 
with results of SIS (solid curve) and NFW (dashed curve) models using the best fit parameters 
in an Q m = 0.3, = 0.7 cosmological model. 
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Fig. 12. — The 68% confidence region in the c-r 2 oo plane for the NFW model prediction to 
the tangential shear of RXJ1347-1145 for Vt m = 1, fi A = (left) and Vt m = 0.3, fi A = 0.7 
(right) cosmological models. The star corresponds to the values that minimize x 2 f° r the 
current weak lensing data. The circle represents the values obtained by Allen et al, 2002, 
from Chandra data. The projections of this error curve onto the c - r 2 oo axes leads to 
estimates of the error in the integrated mass. 
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Fig. 13. — A densitometry measure (5k = R(r < r«) — R(ri < r < r )) of the mass of 
RXJ1347-1145 with NFW (dashed line) and SIS (solid line) models for the best fit parameters 
in the Q m = 0.3, = 0.7 cosmology. Note that the error bars are computed using only the 
number of galaxies in each bin, and are thus non-independent underestimates of the error, 
especially in the central bin. 
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Fig. 14. — The luminosity per unit area in radial annuli in units of solar luminosity per 
pc 2 (triangles) and a best fit scaled mass density that is consistent with gravitational shear 
(solid curve). The dashed curves indicate la variations in the scaled mass density. The 
approximate ratio of mass to light is 91 M /L^ . 



